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Abstract 

We investigate how the qualitative structure of Doppler peaks in the angular power spectrum of the 
cosmic microwave anisotropy is affected by basic assumptions going into theories of structure formation. 
We define the concepts of "coherent" and "incoherent" fluctuations, and also of "active" and "passive" 
fluctuations. In these terms inflationary fluctuations are passive and coherent while topological defects are 
active incoherent fluctuations. Causality and scale invariance are shown to have different implementations 
in theories differing in the above senses. We then extend the formalism of Hu and Sugiyama to treat 
models with cosmic defects. Using this formalism we show that the existence or absence of secondary 
Doppler peaks and the rough placing of the primary peak are very sensitive to the fundamental properties 
defined. We claim therefore that even a rough measurement of the angular power spectrum Ci shape at 
100 < I < 1500 ought to tell us which are the basic ingredients to be used in the right structure formation 
theory. We also apply our formalism to cosmic string theories. These are shown to fall into the class of 
active incoherent theories for which one can robustly predict the absence of secondary Doppler peaks. 
The placing of the cosmic strings' primary peak is more uncertain, but should fall in I » 400 — 600. 



1 Introduction 



The cosmic microwave background (CMB) promises to become one of the most successful bridges between 
theory and experiment in cosmology. As the body of experimental data continues to growjlj], theorists are 
evaluating the impact of this data on the two major paradigms for cosmic structure formation: inflation Q 
and topological defects ||. The so-called Doppler peaks, in particular, have attracted great interest. They 
consist of a system of oscillations, known to be present for most inflationary models, in the CMB angular 
power spectrum C; at 100 < I < 1500. 

The Doppler peaks' height and position have been extensively studied in inflationary scenarios ||, a 
context in which they can be used to fix with some accuracy combinations of cosmological parameters JtJ . 
Although there has been some work on defect models || [J [To| , progress on defect Doppler peak predictions 
has been slow (see however [jll], |l^, |ll|). In two recent letters |l3|, [l9|, we have claimed that regardless of the 
remaining quantitative uncertainties, one could expect dramatic qualitative differences between defect and 
inflationary Doppler peaks. In summary, we found two types of exotic behaviour. Firstly we found that defect 
Doppler peaks may be obtained from the inflationary ones by an additive shift in I. The shift value is roughly 
proportional to the inverse of the defect scaling coherence length. Only very large defects, on the verge of 
violating causality, apply a zero shift. We also found that defects should always soften the secondary Doppler 
oscillations, the more so the larger the additive shift they apply (the smaller their coherence length). Only 
"zero-shift" defects (placing the peaks on the inflationary positions) induce a negligible softening. If the main 
peak is at the isocurvature position one should not obtain more than an undulation. For a smaller coherence 
length (main peak to the right of the isocurvature position) the secondary peaks should be completely erased. 
These two types of phenomena are unheard of in inflationary scenarios. Putting aside "zero-shift" defects, 
we then claimed that measuring the qualitative features of the Ci spectrum at 100 < I < 1500 should decide 
between inflation and defects on fairly general grounds. 

In this paper we explain in more detail the formalism outlined in [[l9) , and present our findings in a more 
systematic form. As in [|l9[ we focus on the basic assumptions of inflationary and defect theories in order to 
isolate the contrasting properties which alone determine the broad Doppler peak features (Section ||) . We 
come up with two concepts: the concept of active and passive perturbations, and the concept of coherent 
and incoherent perturbations. In terms of these inflationary fluctuations are passive and coherent. Topo- 
logical defects are active incoherent fluctuations. Causality and scale invariance are shown to have different 
implementations in theories differing in the above senses. These two concepts then allow us to write down a 
Hu and Sugiyama (H+S) type of formalism 0] tailored for cosmic defects (Sections || and |]). A number of 
simplifications apply, but one also needs to extend the existing work in two ways. Firstly, the way in which 
averages are taken has to be modified for incoherent perturbations. Two approximations may be useful: 
the coherent and the totally incoherent approximations. The factors controlling the applicability of these 
approximations are spelled out in Section 5.4. Secondly the H+S formalism cannot account for radiation 
backreaction, that is, for the radiation itself acting as a source for its own driving force. Causality, in defect 
theories, requires the existence of a large scale nonignorable radiation white noise-spectrum. This leads to a 
backreaction effect which can never be neglected. In Section 3.1.2| we show how to solve this problem. 

We then apply our formalism in two types of study. In Section |5| we perform a qualitative analysis of 
defect Doppler peaks for a generic defect. We explain in more detail the results concerning the position 
of the primary peak and structure of secondary oscillations published in and highlighted above. We 
go beyond the coherent and totally incoherent approximations. We also explicitly translate into C\ spectra 
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arguments initially developed in terms of the photons' energy power spectrum at last scattering. 

On a different front we target cosmic strings (Section^). By using results from a cosmic string simula- 
tion we show that regardless of uncertainties one can safely predict that strings do not exhibit secondary 
oscillations. We also argue that the totally incoherent approximation is applicable to cosmic strings, and 
solve the full H+S algorithm in this approximation. The cosmic string (single) Doppler peak is shown to be 
in the region / « 400 — 600. Our simulation and formalism uncertainties do not yet allow us to make any 
prediction concerning the peak's height relative to the low I plateau. 

2 Review of the Hu and Sugiyama formalism 

We start by highlighting results in j| of which we shall make use. The CMB photons are described by the 
brightness function 0(x, 71,77), defined as the brightness temperature of photons moving in the direction /i 
at point x and time 77. The brightness function satisfies the Boltzmann equation, which may be solved by 
expanding 0(x, /1, 77) in Fourier modes in space and Legendre polynomials in angle. The resulting components 
®i(k,r)) satisfy an infinite hierarchy of ODE's, usually truncated and solved on a computer by a so-called 
Boltzmann code. The angular power spectrum Ci can be obtained from the &i(k, 77) today (77 = 770) by using 

/•OO 

Q = / dfck 2 <|e,(A;,»to)| 2 ) . (1) 
Jo 

where the brackets denote ensemble averages. In the analytical framework developed in [Q, rather than 
making use of a Boltzmann code, evolution is split into three regimes for which analytical solutions are 
found: tight-coupling, recombination, and free-streaming. The idea is that before recombination photons are 
sufficiently tightly coupled to behave like a perfect fluid. For a perfect fluid the infinite series ©z(fc, 77) reduces 
to two non- vanishing components: the monopole Oo (related to the fluid energy density), and the dipole 
(related to its velocity). These satisfy perturbed perfect fluid equations, which can be analytically solved 
using the WKB approximation scheme. After recombination the photons free stream towards us. Although 
the infinite series 0;(fc,?y) is now required to describe the photons an exact solution for its evolution can be 
written. In between the two epochs there is recombination. If recombination were instantaneous one could 
simply glue together the two solutions. However this is never the case. The effects of finite recombination time 
can nevertheless be included in the form of a damping factor D{k), describing Silk damping, an algorithm 
for which is supplied in [0]. A new paper by Battye has recently suggested some potentially important 
corrections to this method ||. The corrections could raise the overall amplitude at small scales, but are 
unlikely to affect the presence or absence of the Doppler peaks, which we address here. 

2.1 Free streaming 

The free-streaming solution is approximately 

©KMo) ~ [Oo + *](fc,77„)j ; (fcA?7») + ©i(fc,77») 2lTl 

rVa 

+ [i> - $}ji(kAr])dri (2) 

where 77, is the conformal time at last scattering, A?? = 770 — 7;, A?/* =770 — 77,, and $ and '5 are the scalar 
gauge-invariant potentials (essentially the Newtonian potential on subhorizon scales). The first two terms 
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represent the free streaming of spatial perturbations (of scale k) at last scattering into angular temperature 
fluctuations (of scale I) here and now. The first term is the free streaming projection of spatial fluctuations 
in the photon energy density (identified as 0o + \&), and the second term the projection of fluctuations in 
the photons' velocity (Doppler effect). The projectors are Bessel functions and peak at I w fcAry*, but with 
some spread. Hence oscillations in the power spectra of 9o + \I r and 81 at last scattering are translated into 
oscillation in the Ci spectrum, but these always appear smoothed in the C; . It is pointed out in 0] that the 
projector of the monopole term is more peaked than the projector of the dipole. Hence oscillations in the 
power spectrum of the monopole Qq at last scattering will be less smoothed out in the C; spectrum than 
oscillations in the dipole Q-y. In addition to the first two terms there is the Integrated Sachs Wolf (ISW) 
term, accounting for perturbations induced on photons in free flight after last scattering by non-conservative 
gravitational potentials. 

2.2 Tight coupling 

The equations ruling the photon fluid during tight coupling are: 

e = ~ei-i 

0r = + + (3) 

where R — § — is the scale factor normalized to 3/4 at photon-baryon equality. In this regime the baryons' 
density contrast At and velocity Vt satisfy the conditions 

A 6 = ^A 7 (4) 

V b = Vy, (5) 

where the photons' density contrast A 7 and velocity Vy are to be found from 

Vy = 01 (6) 

a 7 = i(e + h^j . (7) 

We shall assume here that the main dynamical components are the photons and a defect component (and 
to a minor extent the baryons). Therefore we set the CDM density contrast A c to zero. We also assume 
that Eqn.(^) can be integrated into At = (3/4)A 7 , that is, there are no entropy perturbations. These 
assumptions are not a requirement for the work to be presented, but are merely a simplification. 

It may now happen that the potentials driving the fluid are external, or at any rate, known a priori. If 
this is the case then a WKB solution may be found for the monopole and dipole terms. For large k (say 
k > 0.08h 3 ) this is 

[0o + *](»?.) = [* ~ 4] fa») + * 4 ([cos fcr.fo.) + J(0) sin At. fa.)] [9 + *](0) + Jfo.)) (8) 

©ifa*) = (1 +^)3/4 H 1 + Av*)J(0)][e o + g](0)Binfcr.(fr) 
+ [Jfa) - J(0)][9 + *](0)cosfcr s (? 7>f ) 
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with 



J^lM-^Pd^^G^cosikrM-krM}) (9) 



I{k,ri) = -= / d<n'$(<n')G(ri')sm[kr s (ri)-kr s (ri')} (10) 



G(k, V ) = (l + R)-i/4l-(l + R)l + ^-A (11) 
J(k,7 1 ) = ^ , R (12) 

4fc vT+I 



and where r s (rf) is the sound horizon 



^Tsf^'TlW (13) 



For small fc (say fc < 0.08ft. 3 ) one may simply set R = in the solutions given above. 

The gauge-invariant scalar potentials driving this system may be obtained from Einstein's equations plf 

fc 2 $ = 4Tr(a 2 pA T + p s + 3hv s ) (14) 
$ + * = -8vrfa 2 ^+n s ) (15) 



where a is the scale factor, h = a/a, p (At) is the total matter density (density contrast), and pli/k 2 is 
simply related to the nearly vanishing quadrupole of the photon and neutrino fluctuations. In the tight 
coupling regime one may set II = in a first approximation, and the total density contrast may be written 
as 

pA T = PhA fc + p 7 A 7 = p 7 A 7 (l + R) (16) 
from Equation We have included a defect component with stress-energy tensor 0^ given by 

Ooo = P s 

&ij = p s Sij + {hkj - ^%fc 2 )n s (17) 

as in |22| (but rewritten in terms of Q^ v rather than Equations (|l4|) and ( |l5| ) have the advantage of 

being elliptical, thus allowing a simple connection between sources and potentials. They form a complete 
set as long as the sources are conserved. For the defects this means: 

p s + h('3p s + p s ) + jfcV = (18) 

v s + 2hv s - p s + -k 2 n s = 0. (19) 
3 

3 Topological defects contrasted with inflation 

We now focus on the basic assumptions of inflationary and defect theories and isolate the contrasting prop- 
erties. This will allow us to write down a defect tailored Hu and Sugiyama formalism. We shall use the 
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resulting formalism both in general arguments concerning defect scenarios, and in rigorous calculations for 
particular defects. We define the concepts of active and passive perturbations, and of coherent and incoher- 
ent perturbations. In terms of these concepts inflationary perturbations are passive coherent perturbations. 
Defect perturbations are active perturbations more or less incoherent depending on the defect. 

3.1 Active and passive perturbations, and their different perceptions of causal- 
ity and scaling 

The way in which inflationary and defect perturbations come about is radically different. Inflationary 
fluctuations were produced at a remote epoch, and were driven far outside the Hubble radius by inflation. 
The evolution of these fluctuations is linear (until gravitational collapse becomes non-linear at late times), 
and we call these fluctuations "passive" . Also, because all scales observed today have been in causal contact 
since the onset of inflation, causality does not strongly constrain the fluctuations which result. In contrast, 
defect fluctuations are continuously seeded by defect evolution, which is a non-linear process. We therefore 
say these are "active" perturbations. Also, the constraints imposed by causality on defect formation and 
evolution are much greater than than those placed on inflationary perturbations. 

3.1.1 Active and passive scaling 

The notion of scale invariance has different implications in these two types of theory. For instance, a scale 
invariant gauge-invariant potential $ with dimensions l?l 2 has a power spectrum 

P($) = (|$ k | 2 ) cxfc" 3 

in passive theories (the Harrison- Zeldovich spectrum). This results from the fact that the only variable 
available is fc, and so the only spectrum one can write down which has the right dimensions and does not 
have a scale is the Harrison-Zeldovich spectrum. The situation is different for active theories, since time is 
now a variable. The most general counterpart to the Harrison-Zeldovich spectrum is 

P($) =rfF<s>{krj) (20) 

where Fq, is, to begin with, an arbitrary function of x — krj. All other variables may be written as a product 
of a power of r], ensuring the right dimensions, and an arbitrary function of x. Inspecting all equations it 
can be checked that it is possible to do this consistently for all variables. All equations respect scaling in 
the active sense. 

3.1.2 Causality constraints on active perturbations 

Moreover, active perturbations are constrained by causality, in the form of integral constraints ju], [f5| . 
These consist of energy and momentum conservation laws for fluctuations in an expanding Universe. The 
integral constraints can be used to find the low k behaviour of the perturbations' power spectrum, assuming 
their causal generation and evolution [po[ . Typically it is found that the causal creation and evolution of 
defects requires that their energy p s and scalar velocity v s be white noise at low k, but that the total energy 
fluctuations' power spectrum is required to go like k 4 . To reconcile these two facts one is forced to consider 
the compensation. This is an underdensity in the matter-radiation energy density with a white noise low k 
tail, correlated with the defect network so as to cancel the defects' white-noise tail. 
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We first derive directly the low k behaviour of all variables in our problem. For low k one may set R = 0, 
At ~ A 7 ss 49o, and solve equations (0) coupled to (O) and jl^ ) to find 

e » - ~**irSk (21) 

«-«*W+»«<- - (23, 

where i7 7 = p-f/p, and we have used 3/i 2 = 87173. This solution is in fact valid for all values of k for which 
the dipole ®i is negligible. The lack of superhorizon correlations in the defect network requires the power 
spectrum P(p s ) to have a white noise low k tail. Energy conservation Eqn. ( |l8|) requires that v s also have 
a white noise low k tail. Equations and (^||) show the emergence of the compensation: a necessary 
low k white-noise tail in the radiation power spectrum exactly anticorrelated with the defects' tail. The 
quantity to be cancelled is p s + 3hv s , and not just p s . The density forced to have a k 4 power spectrum 
is U, as shown by (p3|), and not just a combination of Ax and p s . This is precisely the source term for 
the potential which must therefore be white-noise on large scales. Since isotropy requires II s and pll/k 2 
to be constant for small k, the Einstein equations imply that scaling active perturbations produce scaling 
gauge-invariant potentials, which must be white-noise on large scales. In particular P(ty — <f>) = F(krj)i] 3 , 
with F(Q) a non-zero constant. For most realistic defects x A F(x) will then have a single peak, located at a 
value of x = x c corresponding roughly to the "coherence scale" of the defect in question. We will see that 
the place and thickness of the peak in x 4 F(x) are deciding features for the Doppler peaks induced by active 
perturbations. 

These results can be understood by means of integral constraints. Although gauge-invariant in their 
initial formulation, integral constraints have been applied to defects in terms of the gauge dependent energy- 
momentum pseudo-tensor of |l7], [l8|]. Here we go back to their original formulation, and rewrite the defect 
integral constraints in terms of gauge-invariant variables. If a given perturbation has a compact support 
then according to (lj] the following quantity must be conserved 

dV(5T° - x l hST°) . (25) 

It can be checked that the density U = — x l h5T®) is gauge invariant under transformations which go to 
zero outside the perturbation boundary. Since the Universe is unperturbed outside this boundary, setting 
all local perturbation variables to zero in this region is a natural gauge restriction. Assuming a purely scalar 
perturbation, and integrating the second term in ( ^5|) by parts the density hi can then be written as 

U = a 2 pA T + p s + 3hv s (26) 

which is precisely the source term in ( |l4| ) and the quantity given in (|23). Using similar arguments it can 
be shown that the other three integral constraints concern integrals of the form J dV x l U. Using these two 
constraints it was shown in [ pp[ that if the perturbation field is the sum of individual contributions with a 
compact support and which are uncorrelated beyond a certain length, then the power spectrum of the total 
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U field has to be bound by k 4 as k — > 0. This is precisely what we have found by studying directly the low 
k behaviour of the tight-coupled equations. 

The compensation can be implemented in the form of a compensation factor |24j ]. In our language this 
means postulating an equation of state for At of the form 

U = a 2 pA T + p s + 3hv s = lc {p s + 3hv s ) (27) 

Usually one writes 

1 / , 

7r = , 28 

where Xc is the compensation scale. Equation ( |23| ) suggests that if hr) = a then Xc = ^6f2 7 a. Then 
Xc ~ 2.45 in the radiation era. In work in preparation |25| ] we show that j c is more than just a "fudge" 
factor. The radiation backreaction is indeed dominated by a factor of the form of j c . However, secondary 
radiation backreaction effects occur. A dipole driven monopole wave, with k A large scale behaviour and so 
not required by causality, is of particular relevance. It always acts so as to push the compensation scale 
further inside the horizon. Hence we consider a range of Xc going over its limiting value of 2.45. 

3.2 Coherent, incoherent, and totally incoherent perturbations 

Active perturbations may also differ from inflation in the way "chance" comes into the theory. Randomness 
occurs in inflation only when the initial conditions are set up. Time evolution is linear and deterministic, 
and may be found by evolving all variables from an initial value equal to the square root of their initial 
variances. By squaring the result one obtains the variables' variances at any time. Formally this results 
from unequal time correlators of the form 

($(k, V Mk', rf)) = 5(k - k>($(fc, t]))a($(k, r)% (29) 

where a denotes the square root of the power spectrum P. In defect models however, randomness may 
intervene in the time evolution as well as the initial conditions. Although deterministic in principle, the 
defect network evolves as a result of a complicated non-linear process. If there is strong non-linearity, a 
given mode will be "driven" by interactions with the other modes in a way which will force all different-time 
correlators to zero on a time scale characterized by the "coherence time" T c (k, rj). Physically this means that 
one has to perform a new "random" draw after each coherence time in order to construct a defect history 
p3| . The counterpart to ( p9| ) for incoherent perturbations is 

($(k, f?)4(k', r/)) - 5(k - k')P($(fc, rj),rt'-v). (30) 

For \r]' — rj\ = | At/ | > T c (k,rj) we have P($(fc, rj), At7) = 0. For At/ = 0, we recover the power spectrum 
P($(fc, t/), 0) = P{<fr(k,rj)). If there is active scaling we have seen that P($(fc,r/)) = rfF<j>(x). Also the 
correlation 

cor(*( M ),*(M0)= MNt^P^ 

0-($(fc,J7))(7($(fc,T?')) 

must be a function of only x = krj and x' — krj' . Let's call it C$(x,x'). Since there is no time translational 
invariance, this function is not necessarily dependent only on x — x' . Hence, if there is active scaling, we can 
write 

($(k, r,)$(k', rf)) = S(k - k')(r } 3 F*ix)) 1 ' 2 (ri*F*(x')) 1 / 2 C*(x, x'). (31) 
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Detailed understanding of the effects of incoherence on the structure of Doppler peaks requires knowledge 
of the coherence function C(x,x'). However two approximation schemes have been implemented. In one 
effective coherence was assumed so that C(x,x') = 1. This was done in calculations for texture models 
JTT| p"2"| , where coherent statistics were assumed. In another jl3| it was assumed that 

($(k, 77)$(k', rf)) = S(k - k')S( v - V ')P r (<P(k, 77)), (32) 

in which 

P r ($(fc,7j)) = J dAr]P{<P(k,r]),Ar]) (33) 

is the time-integrated power spectrum Q. This is far less naive than just setting C(x, x') = S(x — x') in ( |3l| ) 
and draws on the methods of p4| . It should be a good approximation whenever convolving P($(fc, ??), Arj) 
with functions which vary slowly at the scale of T c (k, r/). Clearly both approximations have drawbacks, and 
it is up to detailed calculations to decide which of the two approximations is better, if any, in each particular 
situation. 

We shall label as coherent, incoherent, and totally incoherent the perturbations satisfying (p9|), (^0|), 
and ([32]) respectively. This feature changes the way the average C l are computed, resulting in a striking 
qualitative difference in the structure of Doppler peaks. 

A comment should be added concerning the possibility of different coherence properties among the 



different components which act as sources for the potentials $ and \P (cf. Eqns. Ql4| ) and (15)). Indeed it 
could happen that say, radiation behaved coherently while defects behaved incoherently. In such a case the 
potentials coherence properties would naturally mix their sources coherence properties in a fashion easily 



predictable from Einstein's equations. However, as we explain in Section 4J, it is a good approximation to 
neglect the a 2 pAx contribution to the potentials except for the compensation. The compensation, on the 



other hand, is given by the large scale solution found in Section |3.1.2| . It should be noted that Eqn. ( [2 2D 
applies for each realization, implying perfect correlation between defects and compensation. It follows 
that the compensation has the same coherence properties as the defect network. Therefore within the 



approximation spelled out further in Section 4.1, the potentials have the same coherence properties as the 
defects. 



4 A defect-tailored Hu and Sugiyama formalism 

In scenarios in which perturbations are driven by a topological defect network, a number of simplifications 
and extensions apply to the Hu and Sugiyama formalism. Making use of active scaling and the low k 
constraints (^lj) to (|24|) one may prove that one may drop all but the convolution terms in (||) and (^|). In 
fact they will be suppressed relative to the convolution terms by a factor of order (riph/v*) 3 ^ 2 > where r] p h is 
the conformal time of the phase transition. Then the multipoles Qi(k,rjo) may be written in the simplified 
form 

ez(Mo) = f 1 'dr ] k<l>(k lV )G(k 1 r,)(D l (k,r 1 ) + V l (k,v)) + 
Jo 

rno 

+ (* -$)(»?*) 31 {kAr}*) + / (*-$)j,(AA»/) (34) 

Jr h 
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with 

T ,l /t x _ 0;-i(fcA^)-q + l)j; +1 (fcA^) 

K(fc,??) _ (2« + l) (1 + ^)3/4 

x [•/(?/*) sin(fcr s (j7*) - kr a (rf)) - cos(fcr s (?7*) - fcr s (r/))] . (35) 

If one assumes coherence then one may simply compute the &i(k, t]q) with this formula, replacing every 
variable by the square root of its power spectrum. By squaring the result one then obtains (\Qi(k, 770 ) | 2 ) , 
providing the C; spectrum by means of ([!]). If the perturbation is incoherent then the H+S formalism must 
be modified at this step. For totally incoherent perturbations (imposing (B2)) one has instead: 



Q = J dkJ\f 1 k 4 P r mk,f 1 )G(k,f 1 ))(D l (k, v ) + V l (k,T 1 )) 2 + 

+ J dkk 2 P{^> - $)(rf)jf(fc A77*) + J dkk 3 [V l ji(kA7 h )a r {<PG)a r («f - $)](»?*) 

dk / dr/k 2 P r (^ - $)jf(kArj) . (36) 

For incoherent perturbations the problem is always more complicated numerically, as one has to perform 
two integrals in time and one in k. 

Even before performing any calculations one may expect striking differences between defects and infla- 
tionary Doppler peaks. For active perturbations the convolution terms dominate over the primordial terms 
in the monopole and dipole terms. Since the phase fcr s (?7*) enters differently in the primordial and convolu- 
tion terms one may expect differences in the Doppler peaks' position. Also, in obtaining the Ci for coherent 
active perturbations one integrates an oscillatory function and then squares the result. Some oscillatory 
structure can always be expected in the C;. For totally incoherent active perturbations, on the contrary, 
one simply integrates the square of an oscillatory function. This contains a DC level as well as an AC level. 
These will compete, leaving it up to other details of the incoherent perturbation to decide whether or not 
there are secondary oscillations. We will expand in great detail on these two points. 

4.1 Possible loopholes 

Eqns. (|l|), (||), (||), (|[), (|l4|), and ( [ili| ) seem to imply that topological defects' C\ spectra may be analytically 
computed just from knowledge of the defect two-point functions 

(0^(77)6^(77')). (37) 

This includes knowledge of time-time correlators, and also cross-correlators involving different stress-energy 
components. However there are two loopholes in this statement which we now spell out. 

In the H+S formalism one assumes that the potentials are external to the photon-baryon fluid, or at 
least, are known a priori. The defect potentials are indeed external, since defects evolve according to their 
own independent dynamics (the so-called stiff defect assumption |l7j]). However we know that we cannot 
neglect A 7 on large scales, as the compensation is of the same order of magnitude as the defect perturbations 
themselves. Since A 7 acts as a source for its own driving potentials we have a loophole in the formalism, 
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as we cannot account for radiation backreaction effects. This is not necessarily a problem, as we have an 
explicit low k solution for the radiation (Eqns. ([H]) to (]24|)). Hence we may account for the nonignorable 
radiation backreaction in the form of a factor 7 C , as in Eqn.(p8|). Thus we include radiation backreaction 
where it is required by causality, and ignore it otherwise. Many of the uncertainties in the backreaction may 
be mimicked by leaving \c as a free parameter. We shall see that the particular Doppler peak features we 
are interested in are insensitive to this uncertainty. In future work [ p3| we will write down solutions for the 
tight-coupling equations driven by a defect source, which incorporate exact radiation backreaction. These 
show that 7 C is in fact the leading order backreaction effect. In another approach J2f| we will simply solve 
the differential equations (||) directly. Although this approach can cope with the compensation exactly, it 
has the disadvantage of requiring a Monte-Carlo technique. 

A second loophole concerns the measurement of the two-point correlators ([37J) from simulations. It may 
happen that we can only measure one term in the overall defect associated stress-energy. For instance a 
full account of cosmic string evolution requires consideration of evolution of tiny loops and the gravitational 
radiation they dissolve into. These are usually overlooked in string simulations. However it is extremely 
important to keep track of every defect stress-energy term, as one needs them all to ensure energy conser- 
vation, without which Eintein's equations are not integrable. If however one assumes that the overlooked 
terms are a minor contribution to the total stress-energy tensor, then they can be neglected in our formal- 
ism. Equations (|l4|) and ( p"5| ) can be solved separately for every stress-energy component. If gravitational 
radiation is negligible then so are its induced potentials <& and 'J. This statement remains true even though 
one needs the gravitational radiation component to integrate the whole set of Einstein's equations. 

These loopholes concern calculations of Ci spectra for concrete motivated defect scenarios (such as cosmic 
strings or textures) . They do not affect the general analysis performed in Section ^|. 

5 Doppler peaks for a generic defect 

We now undertake a qualitative general analysis of peak position and secondary oscillation strength in 
the monopole power spectrum at last scattering. We input qualitative perturbation features, as general 
as possible, for which we output only qualitative peak features. The idea is to find out which novelties 
in defect peak features are generic defect properties as opposed to mere accidents pertaining to particular 
defect models. Also we may learn when and why the coherent and totally incoherent approximations are 
too crude to be useful. 

The free-streaming solution (||) suggests interpreting (0 O + *)(?7*) as the intrinsic anisotropy 5 7 , that is, 
the anisotropy due to photon energy density fluctuations at last scattering. This anisotropy has the sharpest 
free-streaming conversion from k to I and is responsible for the overall Doppler peak features. In particular 
the Doppler peak's position and the structure of secondary oscillations may be inferred from the power 
spectrum P(S- y )(r] Sf ). However it should be borne in mind that oscillations in P(<5 7 ) overestimate oscillations 
in the C/. If we stick to low values of fl^h 2 then the position and structure of the peaks is only mildly 
dependent on these parameters. In this Section we neglect this dependence altogether by implementing an 
approximation where R = 0, r s {rf) — r\j\f%. This is a simplification for illustrative purposes and not a 
formalism requirement. Then, from (||), we have 




(38) 
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Figure 1: The dimensionless power spectrum of the monopole term fc 3 P(Oo + ^(jj,) plotted as a function 
of k* for adiabatic (line) and isocurvature (dashed) passive fluctuations. 

If there is active scaling, in the sense of Eqn.©, one has 

((* - *)(>7)(¥ - $)(?/)) = {r?F(x)f /2 {>q' 3 F{x')) 1/2 C{x, x'). (39) 



with x = kr), x' = krf '. Here F(x) is the structure function (as in |24| ) or scaling factor (as in Eqn.(|2C|)) of 
the potential \& — $, 

P(* - $) = ^Fix) (40) 

and C is its coherence function. The dimensionless power spectrum k 3 P(6y)(r)*) is then only dependent on 
x* = krj* and may be computed from 

fc 3 P(9 + 40(77*) = xlF(x*) - -?=(xlF{x*))V 2 (*' dxC(x, x^)(x 3 F{x))^ 2 sm X * ~ ' 



V3 V * V " Jo V ' " X " V3 
+ ±£*J*'dxdx' (ia^F(x)) 1 ^^^j(x^x , )C(x,^) (41) 

If we use a coherent approximation (C = 1) this becomes 

fc 3 P(9 + mM = ((xlF(x*))V 2 - ^['dx (x 3 F(x))^ 2 sin ^^>j (42) 

If we use a totally incoherent approximation one has instead 

fc 3 P(9 + *)(^) =xlF(x*) + l f dxx 4 F r {x)sm 2 ^=^ (43) 

•J Jo v3 



in which F r is the structure function for the time-integrated power spectrum of \f — $ as defined by ( |32[ ) 
and @: 

P r (*-$) = r) 4 F r (x) (44) 
11 



The rough peak structure in passive theories may be obtained by dropping all but the primordial terms Q , 
thereby ignoring sub-horizon processing. Thus for Harrison- Zeldovich initial conditions one has: 

fc 3 P(9 + *)(?7*) oc cos 2 ^ (45) 

V3 

fc 3 P(9 + #)(J7,) oc sin 2 ^= (46) 

v3 

for adiabatic and isocurvature passive fluctuations. For reference we have plotted these spectra in Figure [j]. 
The maxima are at scales x* m j\rh = run and i^/v3 = (m + 1/2)tt, respectively. There are true zeros in 
the power spectrum, which get smoothed out into valleys in the C; spectrum. These valleys are further 
filled by the dipole, but since this is suppressed, the oscillatory structure survives. The free-streaming 
monopole projector is simply a Bessel function ji(x), which has its main peak at x « I. Then the peaks in 
fc 3 P(0o + ^)(?7*) will be converted into angular scales l m (770 lr\*)x* m . An improvement on this formula 
may be obtained by noticing that the peak of ji(x) is in fact at x ~ I + O.8Z 1 / 3 . 

5.1 A representative generic defect 

We now cut a representative section through the infinite dimensional defect parameter space. The physical 
inputs into ( f4l| ) are the coherence function C(x,x') and the potential structure function. Both may be 
obtained from the sources via Einstein's equation: 

/c 2 ($ - *) = 87r(a 2 pA T + p s + 3hv s + fc 2 IT) (47) 

Then if F p s is the defect energy structure function 

P{p s ) = F 2 pB /r) (48) 

one has 

x*F(x) = (87r) 2 7c 2 F 2 s r 2 (49) 

in which we use a compensation of the form (|27j), and have assumed that the defect's v s and II s are related 
to its energy via equations of state, leading to the factor T. F^T 2 must go to white noise at low k, and 
for realistic defects it falls-off like l/x n at large k. Then F(x) must be white noise at low x, tailing off 
like l/x n+i at high x. x 4 F(x) will then have a peak at a turn-over scale close to x c = 2irrj/^ c , where £ c is 
approximately the coherence length of the defect (the wavelength where F 2 a T 2 switches from white noise to 
power law fall-off). This means that x 3 F(x) and x 4: F r (x) will also have a peak, typically not very far off the 



peak in the potential source structure function x 4 F(x). Inspecting Eqns. ET|, 42, and [43| we see that more 
than the exact form of the potential structure functions, it is the place and thickness of this peak that will 
determine the Doppler peak structure. For the purpose of the qualitative discussion to be carried out we 
then assume for definiteness that for our generic defect 

F(x)=exp^ (50) 

placing the peak of x 4 F(x) at x c = 2a. This form allows us to consider a one parameter family of structure 
functions for which both the peak position and the width are determined by the single parameter a. In the 
familiar cases this "one scale" feature of the structure function is realistic, but it is worth noting that exotic 
cases could provide exceptions to the straightforward analysis which follows from (Eoh. It is an important 



12 



question of principle to know how small x c may be before causality is violated. The value x c w 2.7 was 
suggested in |l5). Recently, a more systematic analysis of this question has been given by Turok|ll| who 
gets similar results. 



Strictly speaking the structure function (50) violates causality. ^From ( \L9\) it implies a structure function 
for p s which once inverted into real space reveals vanishingly small, but not strictly zero superhorizon 
correlations. If one wants to be pedantic about this, one may simply set the real space structure function 
to zero outside the horizon, reinvert to Fourrier space, and find the corresponding modification to F(x). 
Depending on how smoothly the cut-off is done a set of sub-dominant oscillations may or may not appear 
in the modified causal F(x) (as stated in |L6]]). In any case, the correction in F(x) is very small, and it 
certainly does not affect the integrals leading to the C; spectra needed for the sake of our argument. For 
this reason we shall ignore this detail in the rest of our discussion. 

For a compensation satisfying j2^ ) the compensation coherence time is the same as the defect coherence 
time. Hence the potential coherence function is also the defect coherence function. We will assume that our 
generic defect has a coherence function of the form 

CM = expfc^ (51) 

in which the FWHM coherence time is given by 8 C w 2.35r c . 

5.2 The peak position for active perturbations 

If the peak in x i F{x) is sufficiently thin, then from (|4l|), (0), and (|4^), the monopole spectrum peaks will 
be at x* m = (m — l/2)n>/3 + x c for coherent, incoherent, and totally incoherent fluctuations. Then active 
perturbations merely apply a phase shift of value x c — Try/3/2 to an adiabatic type of spectrum (cf. Eqn. fl5|) . 

For x c ~ 2.7 (not impossible, but probably unrealistic because it is very close to the smallest turnover 
point allowed by causality ]l5| , |l6| |) the monopole peaks are at the adiabatic positions. As x c increases from 
the adiabatic position the peaks are shifted to smaller scales. For x c rs 5.4 they are out of phase with the 
adiabatic peaks (as in pi]]). For x c > 8.5 the peaks start only in the adiabatic secondary peaks region. For 
standard values of Q& and h these three cases would place the main "Doppler peak" at I » 230, 350, and 
500, respectively. Therefore the placing of the peaks is not a generic feature of active fluctuations. Active 
perturbations simply add an extra parameter on which the Doppler peaks position is strongly dependent. 
In general we should expect that for the same f2, Qf,, and h, active perturbations will apply to the predicted 
CDM adiabatic peak position a shift of the form 



ttV3 



l^l + ^-f-J (52) 

The secondary peaks' separation is not changed, in a first approximation. This is to be contrasted with 
non-flat inflationary models where Ci(Vt = 1) is taken into Cjq-i/2. The defect shift is additive whereas the 
low-J7 shift is multiplicative, a striking difference that should always allow us to distinguish between low f2 
CDM and = 1 high-x c defects. 

If the peak in x i F(x) is thick then the issue of coherence comes into play but all we have said still 
applies if the perturbation is exactly coherent. This is illustrated in Fig. || using the generic model defined 



in Section 5.1. A formula for the peaks of P(&o + *&){k, ry*) for more general structure functions of coherent 
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Figure 2: Four potential source structure functions (F(x)x where F(x) is the potential structure function) 
and their corresponding spectra fc 3 P(9o + ^)(v*) (dashed lines) assuming coherence. The continuous line 
in the bottom diagram represents a passive adiabatic spectrum. For a sufficiently small £ c one obtains the 
adiabatic structure of peaks. As £ c increases one may obtain the isocurvature peaks. However one cannot 
realistically push the peaks much further to the right without destroying the symmetry of the secondary 
peaks. 

active perturbations can be obtained by computing the structure function "primitives" : 

Fi(ac*) = / dxx a/2 y/F(x) cos(x/\/3) 



F 2 (x*) = / dxx 3/2 ^F(x)sm(x/V3) (53) 



and noting that 



fc 3/2 (6 + <&) = {xlF{x*)) l l 2 - F x {x*) sin ^ + F 2 (x*) cos ^= (54) 

v 3 v 3 



We see that defects mix the adiabatic and isocurvature modes (cf. (^5|) and ([h3])). The defect structure 
function primitives determine the mixing proportions, and therefore the shift applied to the peaks. For peaks 
occurring when x%F(x*) m the peaks' position is given by solutions to the equation 

x* Fi(x*) 

ta,n—= = -—— — - (55) 
V3 F 2 (x.) 
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Figure 3: Four potential source structure functions (F(x)x where F(x) is the potential structure function) 
and their corresponding spectra k 3 P(Qo + V E , )(?7*) (plotted with points) assuming total incoherence. The 
continuous line in the bottom diagram represents a passive adiabatic spectrum. One may obtain (softer) 
secondary oscillations at the adiabatic position for totally incoherent perturbations. As the spectrum shifts 
to the right (larger x c ) the secondary oscillations disappear very quickly. 



5.3 Secondary oscillations for totally incoherent perturbations 

For totally incoherent perturbations with a thick x F(x) peak, secondary oscillations may never show up. 
All we have said on peak position still applies to the main peak position (it is not difficult to rewrite ( [531 ) 
and (^) for totally incoherent perturbations). However secondary oscillations are erased for large x c . The 
situation is illustrated in Fig. || for the generic defect defined in Sec. 5T, using C(x,x') = S(x — x') and 
F r (x) = F(x). Although there are never true zeros in P(Oo + $f)(x*) it is possible to obtain significant 
oscillations if the main peak is at the adiabatic position. In this case the structure function is so narrow 
that the defects only source perturbations for a brief "impluse" time for each scale. However the secondary 
oscillations disappear very quickly as the main peak approaches the isocurvature position, and the structure 
function correspondingly broadens. For x c > 7 there are no significant secondary oscillations. The absence 
of secondary oscillations is therefore not a prediction of incoherence, as even defect's total incoherence is not 
sufficient for the secondary oscillations to be erased. 

The low-x c incoherent oscillations present in P(Oq + ty)(x*) are strong enough to survive in the C;. In 
Fig|I| we show the result of solving the full H+S algorithm for three of the totally incoherent models in 
Sec. 5.1. We have included the monopole, dipole, and monopole-dipole (interference) terms, but dropped 
the ISW term. We took flf, = 0.05, and h — 0.5. For the marginal value x c « 2.7 one may reproduce the 



15 



10 



100 

J 



1000 



Figure 4: The Ci spectrum for three exponential structure functions assuming incoherence and dropping 
the ISW term only. We see that incoherence is not by itself sufficient to erase the secondary oscillations. A 
totally incoherent defect with a primary near the adiabatic position would exhibit secondary oscillations. 



rough sCDM features even with a totally incoherent defect, and with a realistic non-zero value for fit- This 
is somewhat in contradiction with the work of Hu and White p6| , a fact we examine more closely in work in 
preparation ^7j. This example also generalizes Turoks' "confusing defect" |l6|, as it shows that confusing 
defects and sCDM is not only possible, but also that it can be done with defects with any incoherence 
properties. 

These results can be understood from the fact that if the structure function is thin, then each mode is 
only active for a short time. The defect incoherence may not then have a chance to manifest itself in the 
monopole oscillations it drives. This situation may realistically be realized for low x c defects only. If x c is 
high each mode is then active for several expansion times, during which it incoherently kicks the photon 
plasma. The secondary oscillations are then erased or at least softened. 

We may be more quantitative and define the strength of secondary oscillation as the difference between 
power spectra maxima and minima in units of the maximum value: 

_ (k 3 P(6 7 ))(k miV «) - (fc 3 P( ( 5 7 ))(fc n ,^) 

(fe3P(5 7 ))(fe m ,77*) W 

where k m and k n are adjoining maxima and minima in the monopole power spectrum. For passive or coherent 
active perturbations uj m — 1. Approximating the peak by a step function with width 5x c (not too large) 
an estimate of uj m for totally incoherent active perturbations may be obtained. In this approximation the 
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dominant second term in (kt3f) becomes 



and so 



2sin(fa c /V3) 



5x c / \/?> + sin(fe c /v3) 

As 5x c — > one has u> — ► 1, and so the secondary oscillations of totally incoherent perturbations are 
indistinguishable from coherent or passive oscillations. For, say, Sx c rs 5.2 the oscillations are a mere 
undulation with lo = 0.1. 

5.4 Secondary oscillations for incoherent perturbations 

If a perturbation is incoherent, with finite coherence time 9 C , then the strength of its secondary oscillations 
will be overestimated (underestimated) by the coherent (totally incoherent) approximation. In general all 
we have said about peak position applies. The strength of the secondary oscillations will fall somewhere in 
between the coherent and incoherent predictions, the controlling factors being a combination of 5x c and 9 C . 

For a given x c there will be a coherence time 9 C above which the perturbation is effectively coherent, 
and another 9 C below which the totally incoherent approximation is a good approximation. The strength 
of the secondary oscillations w is a good indicator of how good or bad the coherent (predicting u> = 1) and 
totally incoherent (predicting lo as in d58|)) approximations are. We have found the rule that an incoherent 
perturbation is effectively coherent if 8 C ^s> \/25x c , and effectively totally incoherent if 6 C -C \ 28x c . 

The situation for the defect defined in Sec. 5.1 is shown in Figure ^|, where we show spectra in a grid of 



values of x c and 6 C « 2.35r c . In the second and third row we also present results in the totally incoherent 
(dotted line) and coherent (dashed line) approximations. We see that large x c perturbations not only erase 
the secondary oscillations if perfectly incoherent but also require unreasonably large 9 C to deviate from the 
totally incoherent approximation. If x c = 9 (third column of Fig. |E|) then only for 9 C 20 do any meaningful 
secondary oscillations exist. The coherent approximation requires unacceptably large values for 9 C . For all 
reasonable values of 9 C the totally incoherent approximation, on the other hand, simply smoothes out the 
already very soft or nonexistent secondary oscillations. We will see that cosmic strings fall into this category. 

For perturbations with low x c , close to the causality lower bound, smaller values of 9 C are required for 
effective coherence. However these perturbations also show secondary oscillations in the totally incoherent 
approximation, and therefore for all values of 9 C they do not discredit this approximation either. The first 
column in Fig. [5] shows the case x c = 3. For 9 C < 3 the actual oscillations are not far off the oscillations 
present in the totally incoherent approximation. For 9 C ~ 6 one already has lo sa 0.9, and so the coherent 
approximation becomes a very good quantitative approximation. In this regime the real case and both 
approximations do not differ significantly. 

For intermediate values of x c there will normally be meaningful secondary oscillations but peculiarly 
softer than coherent oscillations. Both approximations clearly have shortcomings in this regime. In the 
second column in Fig. |5| we see what happens for x c — 6 (texture type). One needs a rather large coherence 
time (9 C « 20 ) for the oscillations to be as strong as coherent secondary oscillations (say to = 0.9). On 
the other hand the totally incoherent approximation shows nearly non existent secondary oscillations, a 
situation only realized for 9 C < 2. For all the 9 C in between one should do the full calculation to get the right 
qualitative picture. 
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Figure 5: A grid of models as in Sec. 5.1, for various values of x c and 9 C « 2.35r c . The dotted line on the 
second row corresponds to a totally incoherent approximation, and the dashed line in the third row to a 
coherent approximation. 



We may explain these facts semi-quantitatively by repeating for incoherent perturbations the calculation 
leading to the estimate (|58|). The structure function x 3 / 2 F 1 / 2 (x) is then replaced by a step in x centred at 
x c and with width Sx c , and C(x,x') by a step in x — x' centred at zero and with width 8 C . By examining 
the third term in ( |4l| ) one then sees that if 9 C > ^/2Sx c the perturbation is effectively coherent and uj = 1. 
The analogue of ( |57| ) for incoherent perturbations is lengthy and unilluminating. In the limit 9 C <C \/2Sx c it 
may be approximated by the rectangular region of the integration domain: 



6 C /(2V3) I ^ fec/v^ V3 1 



and so 

2 sin (Sx c / \/3) 
Sx c / \/3 + sm(5x c / 



(60) 



IS 



with dx c = Sx c — c /y2. We see that finite coherence time tends to increase the oscillations strength above 
their totally incoherent value, but not by much if 9 C <C V2Sx c . 

5.5 Translation into Cj's 




Figure 6: 1(1 + l)Ci spectra for a grid of models with various values of x c (related to the defect coherence 
length) and 6 C ~ 2.35r c (the defect coherence time). We have included the monopole term (dash) and dipole 
term (point-dash), Silk damping, and free-streaming. The monopole term is always dominant. 

With topological defects, as with inflation, the monopole power spectrum is the dominant term in the 
Doppler peak region. We illustrate this statement in Fig. |6| using a grid of incoherent models as in Fig. [| with 
variable values for x c , and scaling coherence time 9 C ~ 2.35r c . We have included the monopole and dipole 
terms, but dropped the ISW term. We also included Silk damping and assumed that fif, = 0.05, h — 0.5, and 
To = 2.12QK (where To is the CMB average temperature). We see that the dipole (Doppler) term is always 
subdominant, and that the monopole free-streaming further softens the spectrum's oscillations present in 

p(e + *). 

Figure |^ confirms in terms of C; 's what we have inferred from the monopole power spectrum regarding 
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Figure 7: The cosmic strings' energy power spectra as measured from a simulation. This data was used to 
produce the approximate fit to F p *(x) = P{p s {k,rf))rj given in (6f). 

peak position and structure of secondary peaks. The salient features can be summarized as follows. The x c 
parameter controls the peak position. For the extreme value x c « 2.7 the peaks appear on the adiabatic 
position. They suffer an additive shift to the right for larger values of x c . The strength of the secondary 
oscillations depends on both x c (which is connected with structure function width in realistic models) and 
6 C . For x c « 2.7 there are secondary oscillations regardless of the exact 6 C value. This is a confusing defect, 
as not only does it place the Doppler peaks on the adiabatic position, but also the peak structure is quite 
insensitive to the defect incoherence. For larger x c the secondary Doppler peaks survive only if the defect 
coherence time is much larger than x c . This condition seems unphysical for large x c so we expect realistic 
defects with large x c not to have secondary oscillations. 

Here is a rough guide to standard defect theories. Current understanding places the cosmic string models 
on the top right corner of Figure || (large x c , t c smaller than 3). They should have a single peak well after the 
main adiabatic peak. Textures fall somewhere in the middle of the figure (x c around 6, coherence time not 
yet measured). Their main peak should be out of phase with the adiabatic peaks. This is an accident related 
to the x c value for textures, and not a robust defect feature. Texture secondary oscillations should exist 
but be softer than predicted by the coherent approximation (used in jlll [l^]). How much softer depends 
on the exact value of the texture's 6 C . If their coherence time is of the same order as strings (r c ~ 3) their 
secondary oscillation will be very soft. 

6 Do cosmic strings have secondary oscillations? 

Besides the general analysis performed in Section ^ we may also target concrete defect scenarios. The 
parameter space scanning methods developed in Section || may then be useful in allowing use of partial or 
uncertain information obtained from simulations. The idea is to vary the physical inputs of the calculation 
within the simulations' uncertainties, or to fill in what simulations left undetermined in the most general 
form. One may then evaluate the impact of our uncertainties on the final result. It may happen that in 
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Figure 8: The measured cor(//(fc, 40), p s (k, rf)) and the coherence function C p s(x,x') it implies. 

spite of simulation uncertainties some qualitative Doppler peak features are already robust predictions. We 
will argue that this is the case for the absence of secondary oscillations in cosmic string scenarios. 

6.1 What can be measured from simulations 

In Fig. |f]we plot the structure function of the strings' energy density defined by P(p s (k,rj)) — F 2 s (x)/r), as 
discussed in lEBl. There we show that this is well fitted with 



Fl = 



N 



(1 - 2bx c x + bx 2 ) 1 / 2 



(61) 



where N is a constant, b « .006, and x c w 10. Although we cannot measure the i < 1 region with the 
simulation it is known that a white noise behaviour is to be expected. Therefore we may safely take ([H]) as 
a valid extrapolation. The strings' energy coherence function is plotted in Figure. |8[ On the right we plot 



cor(p s (fc,40),p s (fc,7 ? )) = 



a (fi s 40))cr(p s (k, T))) 



(62) 



and on the left the C(x,x') inferred from this measurement. The region x < \,x' < 3 has been left 
undetermined by the simulation, but since C(x 1 x 1 ) is symmetric this leaves undetermined only x < 1, x' < 1. 
We shall fill this gap with the widest range of possibilities. The last input required is T. In the Appendix 
we give some general properties imposed on equations of state by energy conservation. 



6.2 The effects of the compensation and coherence time 

We assume that stress-energy components other than p s can be ignored for the purpose of studying secondary 
oscillations (we take this question up further in the appendix). We therefore set T = 1. The main uncertainty 
from the simulations concerns the strings coherence time. From Figs. ||| this can be estimated to be 8 C w 3 — 6 
(corresponding to r c w 1.5 — 2.5). There is some evidence that r c decreases at large x. We shall fill in the 
region x < l,x' < 1 in three ways: same coherence function as in outer region; incoherent filling, in which 
one sets r c = in this region; and coherent filling, in which one sets C = 1 in this region. These cases 
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Figure 9: The top three rows show a grid of cosmic strings' spectra P(Oo + ^) for r c and Xc varying within 
the allowed uncertainties. In dashed and dotted lines we have plotted a coherent and a totally incoherent 
approximation. The last row addresses the issue of how big the coherence time of a cosmic string would 
have to be for secondary oscillations with ui ~ 0.5 to appear. The value r c = 4 is found, clearly ruled out by 
simulations. 

cover the extreme possibilities for the behaviour of C where it was not measured. Another uncertainty is the 
compensation scale Xc- This may be liberally placed in the region 2 < Xc < 5. Within the framework of all 
that is already known this uncertainty has little impact. 

The top three rows of Fig ^9] show a grid of spectra in which r c and Xc are varied within the allowed 
uncertainties. In dashed and dotted lines we have also plotted a coherent and a totally incoherent approxi- 
mation (using P rather than P r ). We have redone Fig with coherent, incoherent, and trivial fillings in the 
x < l,x' < 1 region and found the results nearly identical. Clearly strings tend to erase secondary oscilla- 
tions. For central values favoured by simulations this is done very effectively. For the marginally acceptable 
case where r c = 2.5, \ = 2 we cannot rule out the existence of very soft undulations in P(6o + '&). The 
bottom row of Fig.^ addresses the issue of how big the coherence time of a cosmic string would have to be 
for secondary oscillations to show up. Even accepting u> — 0.5 as a reasonable oscillation strength, it appears 
that 8 C « 10 is required. This is ruled out by simulations. 

In order to argue that despite simulation uncertainties, we have established cosmic string's absence 
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Figure 10: The marginally acceptable (given simulations) cosmic string parameters which are most favourable 
for secondary oscillations: large coherence time 9 C = 6, coherent filling in the region where C(x, x') could not 
be measured, low values of compensation scale \ c . On the right we show l(l + l)Ci spectra from the monopole 
(dash) and dipole (point-dash) terms, in these cases. As for inflation the monopole dominates. Even pushing 
causality (x c — 2) one does not obtain meaningful oscillations. On the left we show the structure functions 
x A F (x) for the various Xc parameters considered. Broad peaks centred around x = 10 are always obtained. 

of secondary oscillations, we now concentrate on the marginally acceptable scenario most favourable for 
secondary oscillations. We choose the extreme value r c = 2.5 and fill the unmeasured domain of C coherently. 
We then allow Xc to take values 2, • • • , 5. As we have seen the monopole power spectrum then shows an 
undulation. These are not even minima, but mere platforms in the rising spectrum. In Fig. [l^ we show the 
Ci spectrum they translate into. We have included Silk damping, combined the monopole with the dipole, 
and free streamed the result into C;'s. Even in the most extreme case of compensation scale Xc = 2 (pushing 
causality) the resulting Ci undulations are extremely soft. Overall it is clear that the coherent approximation 
for cosmic strings will grossly overestimate the oscillatory structure of the spectrum. The totally incoherent 
approximation, on the other hand, simply exaggerates the lack of oscillations found in the real case. 

This result may be understood by looking at the form of the potential source term x 4 F(x) in Fig. [l0|. 
This is dominated by the shape of F p e(x). A low compensation scale may enhance the power on larger 
scales, and shift the peak in x F(x) to the left. However one would need to push causality limits in order 
to distort the strings' peak significantly, x = 2 is the most that can be done realistically. In general x 4 F(x) 
for cosmic strings will have a very broad peak placed at x = x c « 10. Hence each mode will be active 
for a long time, requiring a very large coherence time for effective coherence. The relatively low coherence 
time measured then suggests effective total incoherence, and since x A F(x) is broad, this will have time to 
erase the secondary oscillations. Also note that the double integral (|4l]) is dominated by a region centred 
at x = x' « 10. Hence it is in this region that knowledge of C(x,x') is most relevant. This explains why 
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our results are so insensitive to the way in which the unmeasured region x < l,x' < 1 is filled. Since we 
also expect that the low x region is where corrections due to string loops and gravity waves are largest, the 
unimportance of this region gives us added confidence in our qualitative results. 

6.3 The cosmic strings Ci spectrum in the totally incoherent approximation 

Since the totally incoherent approximation seems to be justified for cosmic strings we have solved the full H+S 



algorithm in this approximation. We use the expression from 1 24 for the time- integrated power spectrum 



P ^ = TTW (63) 

We consider the two cases (3—1 and (3 = .3 similar to the X and I models in fl24j . We consider only scalar 
contributions. We assume that the defect variables are subject to equations of state of the form p s — r y(x)p s , 
H s — T] 2 j s (x)p s , and v s = r)j v (x)p s . Energy conservation at small x requires that 37(0) = (l/2a) — 1 
and 7u(0) = (1 — 2a)/(3a(4a + 1)), with a — rjh (see appendix). We make use of a string simulation 
to determine the large x behaviour. We find, with large uncertainties, that xj v (x) = s w .1 — .3, and 
x 2 j s (x) = a « .4 — .55. We interpolate between the x ^> 1 and x <§; 1 behaviour. We set II = and assume 
that a 2 pA^ is subdominant except for the compensation and we fix the compensation scale at Xc = 27r. 
Using (|l4|) and jl5| ) we finally obtain the required cosmic strings potential structure functions to be inserted 
in the HS formalism as modified for incoherent perturbations. The results are plotted in Fig. O. The Sachs- 
Wolfe plateau exhibits a "running" tilt ranging from n w 1.4 before I = 10 to n « 1.2 at 30 — 40, although 
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this is quite sensitive to input uncertainties. There is a single Doppler bump located at / « 400 — 600. 
These last two features are remarkably robust against uncertainties in the strings energy structure function, 
compensation and equations of state. The peak position is affected mostly by the energy structure function. 
The lower side of the I ss 400 — 600 range appears to be favoured by X strings, and the upper side by I strings. 
The effect of equations of state on the peak position is negligible. An unnaturally strong compensation on 
small scales could push the peak further to the right but this is unlikely. The absence of secondary oscillation 
is a permanent feature whatever parameters changes one introduces, within the ranges mentioned above. 
The ratio between the peak and the plateau heights, on the other hand, can change by as much as an order 
of magnitude. Even small changes in the equations of state appear to be relevant. This is because the ISW 
plateau and the intrinsic terms probe different combinations of defect stress-energy components. One may 
hope that the peak/plateau height may act as a powerful probe of the defect speed and viscosity composition. 
On the other hand a prediction of the peak/plateau height ratio for strings will have to wait for improved 
simulations. 

7 Conclusions 

We have developed our formalism along two lines. We isolated the differences between inflation and defects 
which we found important for Doppler peak features. These were cast into the concepts of active and passive 
fluctuations, and coherent and incoherent fluctuations. These concepts allowed us to discuss Doppler peak 
features for a generic abstract defect, before addressing any concrete example. In order to address concrete 
examples we had to develop formalism along a second line. We extended the Hu and Sugiyama formalism 
so as to accommodate topological defect theories. The extensions concern mainly the way in which averages 
are taken when the photon fluid is being driven incoherently. We also allowed for photon backreaction effects 
so as to take into account the causality-required compensation. 

We then derived two types of results. We studied Doppler peaks' position and the structure of secondary 
oscillations for a generic defect, and then for cosmic string theories. We found that generic defects place 
the primary peak on or to the right of the adiabatic position (I 220 for sCDM cosmological parameter 
values). The shift to the right, when present, is always additive, in contrast with low Q models, which apply 
a multiplicative type of shift. The value of the shift is controlled by a single parameter x c which can be 
identified as x c = 27r?7/£ c , where £ c is roughly the coherence length of the defect. Very large defects, on the 
verge of violating causality, produce no shift. The smaller the defect coherence length, the larger its shift. 
The texture out-of-phase signature found by jlll [l2| is therefore an accident related to the particular value 
of x c for textures, and not a generic defect feature. 

We also found that the structure of secondary oscillations may be radically different for generic defects. 
This is generally controlled by the ratio of the defect scaling coherence time 9 C and x c . If 8 C x c then 
the defect is effectively coherent, and displays secondary oscillations like the inflationary ones. If this is not 
the case then the secondary oscillations still appear if the defect structure function is sufficiently narrow, 
typically placing the primary peak close to the adiabatic position. Defects with a main peak shifted to 
the right, however, typically have broader structure functions. Provided that 9 C is not larger than x c the 
secondary oscillations get softer, the more so the further to the right the main peak is. Anywhere to the 
right of the isocurvature position, any defect with a realistic 6 C does not show secondary oscillations. 

Given this picture there is good hope for a decisive experiment confronting inflation and defects. Defects 
and inflation can only be confused for f2 = 1 inflation, and a very large defect (x c ~ 2.7). Only this 
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annoying defect would leave no imprint of either its active or of its incoherent nature. Any other defect 
would leave some exotic imprint on the peak structure. An additive shift would ensure no confusion with 
low f2 inflation. Also the secondary peaks would normally appear as softer undulations for not too large x c , 
or totally disappear for any realistic coherent times, if x c is large enough. No inflationary scenario could 
realize these spectra. 

It remains to find out into which region of parameter space (x c ,0 c ) each of the motivated defect scenarios 
fall. This is clearly a quantitative issue to be decided from simulations. We addressed this problem in 
connection with cosmic string theories. We measured the required structure and coherence functions. We 
found that even taking all the uncertainties into account cosmic strings fall into the class of defect theories 
for which the absence of secondary oscillations is a robust prediction. This validates the totally incoherent 
approximation, and we solved the H+S algorithm in this approximation for cosmic strings. The solution 
reveals that given the simulation uncertainties the main peak position should fall in I sa 400 — 600. The 
height of the peak, however, is more sensitive to simulation uncertainties, and so we refrain from commenting 
on it at this stage. 
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Appendix: Active perturbations subject to equations of state 



A complete solution for concrete defect models requires a precise knowledge of the defect stress energy power 
spectra. Knowledge of the stress energy of other related matter components (such as the loops and gravity 
waves into which the cosmic strings decay) is also essential [p8|. Furthermore, the correlations between all 
these components must carefully taken into account. In the case of textures this problem is relatively simple 
because the textures and their primary decay producs are all excitations of the same scalar field which can be 
studied numerically The cosmic string problem is much more complex, and further work is required before 
it is understood clearly. We present here the methodology used for this paper which may, with further 
study, provide a sound basis on which to address these questions. Our starting point is the concept of defect 
equations of state. 

Coherent equations of state 

Active perturbations may be subject to at most two independent coherent equations of state relating p s , p s , 
v s , and II s . Let them be p s = j(kn)p s , and II s = n 2 ^ s (kn)p s for a scaling active perturbation. For a coherent 
perturbation the conservation equations will then impose another equation of state v s = rpf v (krj) p s . The 
equations of state impart perfect equal-time correlation between all four variables (cov(A, B) — <j(A)a(B) 
for any two variables) and so (|2^) becomes the more general statement 

(A(k, vM-k', = (A(-k, r,)B(k', v ')) = 5(k - k')a(A(k, V ))a(B(k' ,r,')) (64) 



where er(.) = y/ P(.). This implies algebraic relations of the form (p s v s ) = <j(p)a(v s ), (pp) = <j(p)a(p), 
or (p 2 ) = (a(p)) 2 . Using them one can square the conservation equations ( |l8| ) and (|l9|), average, take 
the square root of the result, and find that the variances o~(p B ), etc, satisfy the same linear conservations 
equation as the variables p s , etc, themselves. Introducing the variable x = krj, writing a(p s ) — F p s [x)/^ 1 / 2 
and o~(v s ) — F v s{x)ri 1 ^ 2 and denoting d/dx by a prime one has: 

xF' pa + (a(l + 3 7 ) - l/2)F pS + x 2 F vS = (65) 

xfi. + ( 2a + ±) F v * + (\x 2 ls - 7) F p s = (66) 



2/ \3 

where a = rja/a. These can be solved expanding in Taylor series around x = 0. Causality implies that 
F p s (0) = const, and isotropy implies that 7 S (0) = const. To zeroth order equations (|6^) and ( |66| ) imply that 

37(0) = -L-l (67) 

Fixing 7 S (0) can also be done by extending the above calculation up to second order. 

In this context the gauge-invariant formalism has the advantadge over synchronous gauge calculations 
that the potentials are also subject to equations of state of the form $ = n 2r y^,p s and \& = r] 2 j^,p s . This is 
because one may find a complete set of Einstein's equations which are elliptic. One must bear in mind that 
the gauge invariant potentials may be expressed as combinations of synchronous gauge variables, and their 
time derivatives. Hence it is possible to get rid of potential time derivatives in Einstein's equations, assuming 
source conservation, and write equations of state for the potentials. This advantage has one drawback: the 
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defect structure function appearing as a potential source now combines the defect energy with its speed 
and viscosity, whereas in the synchronous gauge approach one only needs its pressure. Depending on the 



information available on the defect this may or may not be a problem. ^From Einstein's equations (14) and 
( |l5| ) we have 

7* = ^i+J^ (69) 
7* = -4-(7ci^ + 27,). (70) 
These allow us to write the required relations 

„(* G ) = ,3/ W ,= ^ J J^ m ( 1 - {1 + R) ^ + §-A (7!) 



Using the conservation equations ( |18| ) and (|19j) we could also write $ and ^ in terms of defect variables (and 
not their derivatives) and find a similar equation of state for the potentials appearing in the ISW term. The 
result is lengthy but straightforward to derive. 

Small scale equations of state 

On small scales the defects cannot be subject to coherent equations of state. For realistic sources it happens 
that Fp3(x) oc x~ n , and 7 and k 2 H s / p s = x 2 j s tend to a constant for i» 1. One may check that equations 
( |65| ) and (|6^) would then lead to over-restrictive constraints (like n = 2a — 3/2). We find from simulations 
that rather than relations like (|64|) on small scales one has 

cor(A(k, 77), B(-k\ 7/)) « (72) 

where A and B are any two defect stress-energy components. This means that the different defect stress- 
energy components are in fact uncorrelated random variables, for which it makes more sense to postulate 
equations of state of the form: 

(P s2 ) = 1 2 (P S2 ) 
(v s2 ) = V 2 1 2 (P S2 ) 

<n s2 ) = v 4 i 2 s(p s2 ) (73) 

We have found for cosmic strings that for large x 

x^j v {x) = 8W.1-.3 
x 2 1 s {x) = crw.4-.55 (74) 

The effect of the equations of state on cosmic strings coherence properties 

^From what we have said above the low x behaviour of the T factor (defined in ([l9])) is 

1 - 2a ' 2 



+ ^+7,(0)) (75) 
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whereas its high x behaviour is of the form 



T z = l + 9a z —+a 2 (76) 

x z 

In both regimes the T factor is of order 1. Therefore T may never modify the shape of the structure function 
x 4 F(x) as given by (|4^), which is determined by F p s and the compensation factor. For this reason we simply 
set r = 1 in our discussion of strings coherence features. On the other hand no doubt T will affect the 
relative normalization of the structure functions used in the intrinsic and the ISW terms, which use different 
equations of state. We have found the the relative height of the peak and low I plateau reflect sensitively 
the defect speed and viscosity composition. 

Incoherent perturbations 

If the perturbation is not coherent but is still subject to two equations of state, then all we have said still 
holds for the power spectra of the perturbation, but not for its time-integrated power spectra. One may try 
to relate power spectra and time-integrated power spectra. Suppose that one can factorize the unequal-time 
correlator as 

P(p s (k, v ),A7 ] )^<j(p s (k,7 1 )) f j(p s (k,r 1 + A7 1 ))C(k,A7 1 ) (77) 

for, say, p s . C(k,Arj) is a coherence function not necessarily symmetric about zero in A?/. Then wherever 
C ~ 1 there is coherence, which is then cut off as C goes to zero for Arj > t c . Suppose that this transition 
is very abrupt (e.g. C is nearly a step function). Then we can expand the time-integrated power spectrum 
written as in (t53]) in the form 



r 0) 

P r (p s (k, V )) = r^Pip^k^+r^aip^k^Wik^)) + J^a( p »(k,ri))a(p s (k,r,)) + 0(Tt) (78) 
where 

/oo 
dAri C{k, A^A^y- 1 . (79) 
— oo 

We have r™ = C(r™). However, if C is symmetric about zero all the even n vanish. We shall assume that this 
is not the case, so that the term in is the dominant correction to order 0{t 2 ). Using the conservation 
equations we can then write 

P r (p s (k, V )) = r c P(p s (fc,r ? ))(l - fP(a(l + 3 7 ) + x 2 lv )) (80) 

^ (2) (2) 

with Tc = Tc /(t c t]). For high k the behaviour of P and P r may be very different. For instance, for cosmic 
strings it is known that P r (p s (k,r))) oc k~ 2 whereas P(p s (k,ri)) oc fc _1 (J2J]). However for small k equation 
(pOf) tells us that that P r and P are related simply by a multiplicative constant. It can be checked that this 



is true for the power spectra of any variable. Moreover the multiplicative constant is the same (r c (1) ) for all 
variables if one ignores corrections in f c . Therefore to this level of approximation the large scale behaviour 
of P r can be simply inferred from ( |67| ) and (|6^). If, however, one is to keep corrections of order fi 2 ^ then 
the situation is more complicated, as the multiplicative constant depends on the variable. For a quantity 
related to p s by an equation of state of the form T(kr])p s one then has for x — > 

P r {T(k V )p s (k, V )) - r c T 2 (k V )P(p s (k,ri))(l - f( 2 )(a(I + 3 7 ) + r'/r)) . (81) 
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~ (2) 

Higher derivatives of the 7's at zero would be required to C(f c ). How important these corrections are 
remains to be determined. 

Modelling the P r on small scales is again more complicated and defect dependent. If the 7's are always 
slowly varying one may use an approximation where formally P r = rjT^a r2 , where = T c /rj, and <j r means 
the value of a as computed for coherent fluctuations but with P(p s ) replaced by P r (p s ). Then the required 
potential properties can be found from 

p r ($G) = r c y 7 L^ 2 (82) 

P r (*-$) = r^l w F r p ? (83) 
with ji ri t and ji SW are the equations of state for the potentials appearing in the intrinsic and ISW terms. 



We have used this approximation in the calculation in Section 6.3 of the Ci spectrum of cosmic strings in 
the totally incoherent approximation. 
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